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Abstract 

With the aid of the vector Ward-Takahashi identity we derive a general formula for the quark 
number susceptibility (QNS) which expresses the QNS as an integral expression only involving the 
full quark propagator at finite temperature and chemical potential. The QNS at finite temper- 
ature and zero chemical potential is calculated with the dressed quark propagator in the Hard- 
Thermal-Loop (HTL) approximation. A comparison of our result with the results of QNS in HTL 
approximation in previous literatures is given. 
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Analysis of fluctuations is a powerful method for characterizing the thermodynamical 
properties of a system. It is believed that enhanced fluctuations are an essential characteristic 
of phase transitions measure of the intrinsic statistical fluctuations in a system close 

to thermal equilibrium is provided by the corresponding susceptibilities. Among various 
QCD susceptibilities, which are associated with the response of the QCD vacuum to external 
sources, the quark number susceptibility (QNS) plays an important role in identifying the 
critical end point in the QCD phase diagram {2, 7-12]. To date, there has been a great 



amount of related works for calculating the QNS, for example, lattice QCD simulation 

■"n nn 

16| . the Nambu-Jona-Lasinio model |5|, 



13 



17 



181 ] and Polyakov-Nambu-Jona-Lasinio 



|20| . the QCD Dyson- Schwinger equation approach 



12 



21 



23| , the hard- 



model (PNJL) 

thermal- loop (HTL) approximation |24j-|27|. and so on. In this paper we try to revisit the 
calculation of the QNS in the HTL approximation by a direct method. 



The partition function of QCD at zero temperature and finite chemical potential reads 



W = JVMDA op{-S[J, 8 M]+/d*x/«(x)7«W}, 



(1) 



where S[q, q, A] is the standard Euclidean QCD action with q being the quark field with 
two flavors (we confine ourselves to the two-flavor case with exact isospin symmetry and set 
A*u — A*d — A 4 ) where \x u and are the chemical potential of the up and down quarks) and 
three colors. Here we leave the ghost field term and its integration measure to be understood. 
The pressure density Vip) is given by 

1 



V(ji) = - lnZ[/i], 



(2) 



where V is the four-volume normalizing factor. From this one immediately obtains the 
quark-number density 

&P(jj.) 1 1 dZ[/j] 



dfj, V Z\p] dfj, 
1 / VqDqVA j d A xq(x) r y4 : q(x) exp {—S[q, q, A; fj]} 



V j VqDqVA exp {-S[q,q, A; fi]} 

where S[q, q, A; fj] = S[q, q, A] — J d A x liq{x)^i^q{x). 



(3) 



On the other hand, the full quark propagator at finite chemical potential is defined as 

j VqDqVA q^xjq^y) exp {S[q, q, A; fi]} 



Gij[fj,}(x,y) 



J VqDqVA exp {—S[q, q, A; //])} 



(4) 
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From Eq. (4) it is easy to obtain the following 

JVqVqVAJd 4 xq(x)-f A q(x)exp{-S[q,q,A;fj]} 
* {GM74} = HP e*p{-S[J,„4ri} • (5) 

where the notation Tr denotes trace over the color, flavor, Dirac and coordinate space indices. 
Comparing Eq. (3) with (5), we obtain a well-known result 



M = -yTr {GM 74 } = -N e N f J ^tr {G^Kpfr*} , (6) 

where N c and Nf denote the number of colors and of flavors, respectively, and the trace 
operation tr is over Dirac indices. From Eq. (6) it can be seen that the quark-number 
density p(fi) is totally determined by the full quark propagator at finite chemical potential. 
It should be noted that without loss of generality, in deriving Eq. (6) we have used the 
unrenormalized language. For the renormalized formulation please see Ref. {29 ]. 

The QNS is defined as the derivative of quark number density with respect to the quark 
chemical potential 

Here, in order to show the difference between our approach and that in the previous litera- 
tures, let us briefly recall the approach in the previous literatures. Substituting Eq. (6) into 
Eq. (7) and adopting the identity 

dG{p,fi) dG-\p,fi) 

— 0^ — = (-)G(p, fj) ^ G(p, fj), (8) 

we have 

X = N c N f J ^tr [g(p, n) 9G ~^ /i) G(p, „) 74 ] . (9) 
Recall that the well-known Ward identity 

ir>,o) = ^M, (io) 

where p denotes the relative momentum of the vector vertex and the corresponding to- 
tal momentum vanishes. Note that at finite chemical potential, the fourth component of 
momentum p 4 = p 4 + ifi. Then 

(-)r l(p ,0;,) = ^^. (11) 
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Putting this equation into Eq. (9) and replacing the integration over the fourth component 
of momentum with explicit summation over Matsubara frequencies, we get the QNS at finite 
temperature and chemical potential in the imaginary-time formalism 

+ °° r d 3 n 

X (T,fj) = (-)N e N f T £ /_^tr[G(p)r4(p,0)G(p) 7 4], (12) 

where p = (p,u n + ifi) with fermion Matsubara frequencies co n = (2n + 1)ttT. Thus a 
model-independent closed integral formula is obtained, which expresses the QNS in terms 
of the full quark propagator and the vector vertex, both of the latter objects being basic 
quantities in quantum field theory. From Eq. (12) it can be seen that in order to calculate 
the QNS by means of this formula, one needs to know both the full quark propagator and 
the vector vertex at finite T and \x in advance. For example, in previous literatures 24j-|27|, 
when those authors calculate QNS in HTL approximation by means of Eq. (12), they have 
to calculate the quark propagator and the quark-meson vertex in the vector meson channel 
at zero total momentum separately. And they have to further check whether the obtained 
quark-meson vertex and quark propagator in HTL approximation satisfy the vector Ward- 
Takahashi identity (as was pointed out in Ref . |27| , in a self- consistent calculation of QNS in 
HTL approximation it is very important to ensure that the obtained quark-meson vertex and 
quark propagator in HTL approximation satisfy the vector Ward-Takahashi identity). In 
fact, one can use the vector Ward-Takahashi identity to obtain the full vector quark- meson 
vertex at zero total momentum directly from the full quark propagator (see, for example, 
Ref. |s0|), so it is not necessary to calculate the vector quark- meson vertex separately and 
check the validity of the Ward-Takahashi identity. As will be shown below, with the aid of 
the Ward-Takahashi identity it is not difficult to show that in the calculation of QNS one 
only needs to know the full quark propagator at finite T and \i. 

Substituting Eq. (6) into Eq. (7) and replacing the integration over the fourth component 
of momentum with explicit summation over Matsubara frequencies, one can obtain a general 
formula for the QNS at finite T and \i 

+00 . j3- a 

X (li,T) = —N c NfT £ J ^— tT [G[»,T}(p,u;M. (13) 

Here we note that at finite temperature and chemical potential, the fourth component of 
momentum pi = uj n + Making use of this fact, Eq. (13) can be rewritten as 

xfaT) = —iN c NfT £ y ^_ta[GM(S WB ) 7 4]. (14) 



From this it can be seen that the QNS at finite chemical potential and temperature is totally 
determined by the full quark propagator G[fi,T]{p,oj n ). At this point some discussions are 
necessary. It should be noted that in the CJT (Cornwall- Jackiw-Tomboulis) or $ functional 
approach, the grand canonical potential is not given by the propagator alone: an infinite se- 
ries of irreducible diagrams is additionally needed. However, in some cases, the situation will 
be changed. For example, it can be shown that the nontrivial \x dependence of the partition 
function of QCD at finite chemical potential and zero temperature is totally determined by 
the full quark propagator at finite chemical potential and zero temperature (more details 



can be found in Ref. 



29(|) . As far as the present work, we do not claim that the partition 



function of QCD at both finite chemical potential and finite temperature is determined by 
the quark propagator alone. We just try to propose a direct and model-independent formula 
(Eqs. (13,14)) for calculating the QNS at finite chemical potential and temperature from 
the full quark propagator. In principle, once the exact form of G[fi,T](p,uj n ) is known, one 
can have a thorough understanding of the QNS at finite chemical potential and tempera- 
ture. However, at present it is very difficult to calculate G[fi,T](p,u n ) from first principles 
of QCD. So when one uses formula (13) or (14) to calculate x(^jT), one has to resort to 
various nonperturbative QCD models. In this paper we shall use the HTL approximation to 
calculate the QNS at finite temperature and zero chemical potential by means of the above 
proposed direct method. 

It is generally believed that when the temperature T is high enough, the HTL approxima- 



tion is a good approximation for QCD [3lM34|. Based on this idea, the authors of Ref. [24] 



first calculated QNS in HTL approximation. Later the calculation of QNS in HTL approxi- 
mation was further studied in Refs. 25M27I] . Here we would like to stress that the purpose of 



our work is to give a direct method to calculate the QNS. In this method the QNS is totally 
determined by the full quark propagator at finite T and //. And in our work we choose the 
HTL quark propagator as an example to illustrate how to apply this method to calculate 
the QNS and compare the obtained result with the existing results of QNS in the literature. 
Here we also would like to stress that in the calculation of QNS in HTL approximation done 
in our work the vector Ward-Takahashi identity is automatically satisfied. 
In the HTL approximation, the quark propagator reads 

GM = -DM— D.(p) 2 ■ < 15 > 



5 



where p = (p,P4), P = p/\p\, P4 = (2n + l)irT (n G Z) are the fermion Matsubara frequencies 
and 

,2 r 



TYl 

D±(p) = -%Pa ± \p\ + -j^ 



ip 4 



ip 4 



Qo[tt\tQi[t^ 



\p\ 



\p\ 



(16) 



In Eq. ffl6|) . m q = gT/\/6 is the quark "thermal mass" with g being the strong coupling 
constant. Qq and Qi are Legendre functions of the second kind. Just as was shown above, 
we have two equivalent formulas (Eqs. (13) and (14)) for calculating the QNS at finite 
chemical potential and temperature. In actual calculation of the QNS using which formula 
depends on which formula is more convenient. If one tries to use Eq. (13) to calculate the 
QNS at finite temperature and zero chemical potential in HTL approximation, one should 
use /^-dependent HTL quark propagator (for very small values of fi) and then take the limit 
fi — > after d/dfi operation in Eq. (13). It is generally believed that when T or fi is high 
enough, the HTL or Hard-Dense-Loop (HDL) quark propagator is a good approximation 
for QCD. So far we do not have a reliable quark propagator at high T and small /i. Hence 
it is not convenient to calculate the QNS at finite temperature and zero chemical potential 
in HTL approximation using Eq. (13). However, in Eq. (14) the d/dfi operation has been 
replaced by d/dpi operation. So, if one tries to calculate the QNS at zero fi and finite T in 
the HTL approximation, one can directly set fi = in Eq. (14) and then input the HTL 
quark propagator at fi = (Eq. (15)) into Eq. (14). By this method one obtains a formula 
for the QNS at zero fi and finite T in the HTL approximation: 

d 3 p (dG 



= ~2tN c N f TY: J 



(2tt) j 
d 3 p 
(2^)3 



dp a 
D' 



2iN c N f T 



drp 
(2^)3 



£4 



(17) 



where ' means d/dp^. From Eq. f )16p we can easily obtain the following 



D' 



D' 



-i + 



\P\ 



rnq 
\P\ 



1 " it) Q '°~ m Q \ 



\p\J 



\p\ 



(18) 
(19) 



To calculate the frequency sum in Eq. ( fTTj) . let us consider the following contour integral 



\mp4 



-3ttT 



--/ 



c 2 



3ttT 
— • — 



5ttT 



FIG. 1: The integral contours in complex p, plane 



where the integral contours C\ and Ci are shown in Fig. 1. The contours are carefully 
selected because (D' + /D+ + D'_/ D 2 _) has a branch cut from —i\p\ to i\p\. When g — y 0, the 



contour turns into the usually emp 
field theory (see, for example, Ref. 



oyed one which can be found in textbooks of thermal 



35|) 



Applying Cauchy's theorem one has the following 



dp, fD\. 

dUC 2 27T \Dl 



D'_ 



— tanh ( 
2 V2T 



Di 



+ V- tanh 
Dl D 2 _ ) 2 V2T 



(21) 



'-i|p| 27T 

To evaluate the integrals on the two sides of Eq. (21), we first use integration by parts to 
obtain 



dp* (&+ 
duc 2 2ir \D% 



- tanh ('-^-\ 
2 \2T ) 



dp, 

C1UC2 27T 



exp(zp 4 /T) 



D_) T[exp(i Pi /T) + lp 



(22) 



and 



M dp, 

i\p\ 2tt 
Ip1+o+ 



Disc 



+ ^ V tanh 
Di D 2 j 2 V2T 



-i|p|+0+ 



dp 4 
27 



D' 



D 2 



— r I - tanh ( 
D 2 /2 V2T 



i|p|+0" 



dp 4 

27 



-i|pl+0" 



Dl 



— 7T j - tanh 
D 2 /2 V2T 



(23) 



7 



j|p|+0+ 



dp 4 ( 1 1 \ i exp(ip 4 /T) 
27 \D^ + iT) T [exp{m/T) + l] 2 



i|p|+0" 



dp 4 / 1 1 \ « exp(ip 4 /T) 
27 ^ + IT J T[exp(^ 4 /T) + 1] 2 ' 



-i|p|+0+ v -i|p|+0" 

Here we note that when Eq. (23) is obtained by integration by parts, the contributions from 
end points cancel each other. 

In the complex p 4 -plane, the function 



exp(ip 4 /T) 



(24) 



V D+ D-J [exp{ip 4 /T) + l] 2 
has four poles and a branch cut from —i\p\ to i\p\ in the imaginary p 4 axis (See Fig. 1). 
These four poles are located at z k = ±iou±(\p\), k — 1, . . . ,4 and w±(|p|) (u;±(|p|) > \p\) are 
determined by the following equations 



\p\(u+ - \p\) 



mi 



|p](w_ + \p\) 



- 1 



+ 1 



l(\-^±\ In U+ + |P1 



m„ 



1 + W ~ | In °°~ + 



(25) 
(26) 



By closing the integral contour Ci and C 2 by large half-circles (note that the statistical 
factor is bounded on the circle) we can get the following 



if dp 4 

7f / IT- 9 {Pi) = 

T Jciuc 2 2tt 



2ixT 



-2ni Res[g(z k )} \ . 
k=i ) 



Using 



Res 
Res 



1 



D+(-iw+) 
1 



= Res 
Res 



1 



.w 



1 



2m 2 



.uj_ 



the sum of residues in Eq. (27) is calculated to be 



2m 



S! = iz_(ifl), 



h ^ " l [exp(c + /T) + I] 2 + [exp(u;_/T) + 1]= 



Therefore we obtain 

i f dp^ 



Ic 



2i j Z + exp(u + /T) Z_ exp(w_/T) 



T J c*iuc*2 2vr * T \ [exp(w + /T) + l] 2 [exp(w_/T) + l] 2 

The integral in Eq. (23) can be calculated to be 



(27) 

(28) 
(29) 

(30) 
(31) 



Disc 



IpI dp A 
2n 

\p\ duo 



exp(ip 4 /T) 



1 1 \ % 
D^ + lT)f [eMm/T) + l] 2 



% fW\ du, . . ... exp(o;/T) 
-f U 2T [P+(W) + P - ( " )] [exp(c/T) + l] 2 ' 



(32) 



where p±(u, \p\) are the familiar spectral functions of the quark propagator (x = cj/\p\) 



P± (u,\p\) = 27r[Z ± (\p\)5(u;-u ± (\p\)) + Z^\p\)d(u + u^\p\))} 
Tim 2 (\ =F x)9(l - x 2 ) 



+ 



2 4 

7T m„ 



m„ m„ , 1 | x 

l T i + ^±^(l T x) In 



1 — x 



(33) 



Combining Eqs. (l2Tj) . ( )3T|) . (|32|) . we obtain the following 





D'_\ 


2i f 






~ T { 



+ 



i r\p\ dui r , . . ex 
+ t/- Ip1 ^ I ' ,+ M + ' , - (w)1 1^( 



exp(w/T) 
w/T) + 1]' 



2i 



where np is the Fermi distribution function 



— [Z + n F (u + )(l - n F (u+)) + Z_n F (uJ)(l - n F {uJ))\ 
\p+(u) + p-(u)]n F (u)(l - n F (u)), 



T 

+ 



i f\p\ du 



(34) 



n F (u) 



exp(u/T) + 1 

Substituting Eq. (134)) into Eq. (|T7|) . we obtain the QNS in HTL approximation 



(35) 



X (T) = Xq {T) + X l{T) 
AN c N f f d 3 p 



2m\ 



T J (2vr) 3 
2N c N f f d 3 p f\P\ du 



u 2 ~^ oj 2 ~^ 

" + ' n F (u+)(l - n F (u + )) + ~ n F (u-)(l - n F (u-)) 



2m 2 



[P+H + P-(u)]n F (u)(l - n F (u)), 



T J (2vr) 3 J-\p\ 2vr 
where Xq(T) is the pole-pole contributions from the quasi-particle mode 

d 3 p 



(36) 



X q (T) 



AN c Nf 



(2tt)' 



n F (u + )(l - n F (u+)) + ~ - n F (w_)) 



2m^ 



2m^ 



and Xl(T) is the cut contributions from the Landau damping 

2N C N } r d 3 p [\p\ du 



Xl(T) 



[p+(u) + p^{u)]n F (u)(l - n F (u)). 



(37) 



(38) 



T J (2tt) 3 J-\p\ 2tt 1 

It is obvious that the obtained QNS contains the contributions from both the Landau damp- 
ing and the quasi-particle poles. This result is quite different from that of Ref. 25]. In Ref. 



9 



251 ] the authors argue that "the cut contributions due to space-like quark momenta do not 



contribute because of the number conservation" . In our opinion, although the number con- 
servation does hold, the quasi-particle mode and the Landau damping mode are excited at 
different energy levels and therefore have different statistical factors. As a result, the final 
result contains both two contributions. Furthermore, even if one ignores the fact that the 
cut contribution from the Landau damping is nonzero, our result for the contribution from 
quasi-particle mode Xq(T) still disagrees with the result of Ref. |25| (please compare Eq. 
(35) of Ref. 25) with Eq. (37) in the present work). Here we should note that various 
approximations employed in the literature lead to different results (see. e.g., Ref. 26)). 
which shows the sensitivity of the quantity under consideration and points to the need of a 
concise approach to arrive at a unique result. This is one of our motivations of doing this 
work. 

The strong coupling constant g(T) can be expressed as 



g(T) 



4:71 X 127T 



(39) 



V (33 - 2N f ) ln(Q 2 /A§) ' 
where A = 200 ~ 300MeV and Q is the momentum scale. In numerical calculations, we 
take Nf = 2, Q = ItxT 36( and the phase transition point T c = 0.49Ao (about 98 ~ 147 
MeV). 
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FIG. 2: x/Xfree as a function of T/T c , where Xfree = N c NfT 2 /3 is the QNS of the free quark gas 



The numerical result is shown in Fig. 2. Just as is shown by Eqs. (36)-(38), the QNS 
has contributions from both quasi-particle mode and Landau damping. In order to show 

10 



clearly these two contributions, in Fig. 2 we plot the contribution from the quasi-particle 
mode and the contribution from Landau damping separately. From Fig. 2 it can be seen 
;hat if one only considers the contribution from quasi-particle mode, as was done in Ref. 



25] , then the QNS calculated in HTL approximation will be smaller than the corresponding 



susceptibility of free quark gas in the whole range of T, and when T — > +oo the quasi-particle 
part will tend to the free quark gas result which can be easily seen from Eq. (37) (when 
T —7- +oo, Z + — y 1, Z-. — > and uj± — > \p\). Fig. 2 also shows that the contribution from 
Landau damping is important even for T ~ lOTc. For the free quark gas, there are only 
particle excitations, whereas in the calculation of QNS in the HTL approximation done in 
this paper there are both quasi-particle excitations and Landau damping mode excitations. 
It is just the fact that the contribution of Landau damping is important which renders the 
QNS calculated in HTL approximation larger than the QNS of the free quark gas. 

Here we should notice that the result in Fig. 2 is different from the lattice result in 
It also differs from some other existing results within the HTL approximation 




and in the PNJL model 



19 



201 ] . The origin of the difference between the result of 



QNS in Fig. 2 and those in the previous literatures can be understood as follows. From 
Eq. (17) it can be seen that in order to calculate the QNS in HTL approximation, one 
needs to know the HTL quark propagator in the whole momentum range. It is well-known 
that the HTL approximation is only valid for external momentum much smaller than T. If 
one ignores this fact and assumes that the HTL quark propagator (15) is applicable in the 
whole momentum r ang e (this is in the same spirit as the one for calculating QNS in HTL 



approximation 25|, |27|]), then one will obtain the QNS result shown in Fig. 2. Therefore 



the result in Fig. 2 is a consequence of ignoring the the range of applicability of the HTL 
quark propagator (15). So in a consistent calculation of QNS in the HTL approximation by 
means of Eq. (17) one should introduce a momentum cutoff A^tl below which the HTL 
quark propagator (15) is applicable. In our calculation we choose Ahtl = gT. We make 
this choice because gT is an important energy scale to identify the soft momentum in HTL 
approximation. If we limit the range of integration of the contribution of Landau damping 
in Eq. (38) to the region |p| < gT, we will get the numerical result of the QNS shown in 
Fig. 3. From Fig. 3 it can be seen that the QNS under HTL approximation is smaller than 
the QNS of free quark gas in the whole range of T and the contribution of the quasi-particle 
mode is much larger than that of Landau damping. Here we should point out that the result 



11 



in Fig. 3 is in fact not sensitive to the variation of the cut-off gT. For example if one set 
the cut-off to be 2gT one would find the change of \ is about 10%. This is because setting 
cut-off in the order of gT would result in xl « Xqi an d consequently, the quark number 
susceptibility is determined essentially within the given approximations by the quark quasi- 
particle contribution. Therefore one would find that the result of x{T) shown in Fig. 3 is 



in agreement with the results in Refs. 



16 



341 ] on the level of 10 percents. In addition, it 



should be noted that our results are gauge independent due to gauge invariance of the HTL 
quark propagator. We also want to point out that in the calculation of QNS in the present 
paper we have used HTL quark propagator at one-loop order. In fact, the method employed 
in the present paper can be easily applied to the calculation of the QNS using HTL quark 
propagator up to two-loop or higher order. 
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FIG. 3: x/Xfree with range of integration of contribution of Landau damping limited in the region 
p < gT 



To summarize, in the present work we derive a new and model-independent formula 
(Eq. (14)) for calculating QNS, which expresses the QNS totally in terms of the full quark 
propagator at finite temperature and chemical potential. Then we choose the HTL quark 
propagator at zero chemical potential (Eq. (15)) as an approximation to the full quark prop- 
agator to calculate the QNS at finite temperature and zero chemical potential by means of 
our method (see Eq. (17) in our paper). As is shown in Eq. (17), when one calculates the 
QNS at finite T and zero fi, one should integrate the momentum from zero to infinity. How- 



12 



ever, it is well-known that HTL quark propagator is only applicable for external momentum 
much smaller than the temperature T and one should choose a boundary value Ahtl below 
which the HTL quark propagator is applicable. In this paper we have chosen A HTL = gT 
because gT is an important energy scale to identify the soft momentum in HTL approxi- 
mation. When limiting the range of integration of the contribution of Landau damping to 
be p < gT, we find that the result of QNS calculated by our method is consistent with the 
results of previous literatures. 
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